A Comparison of Gene Expression Profiles between Glucocorticoid Responder and Non-Responder Bovine Trabecular Meshwork Cells Using RNA Sequencing

The most common ocular side effect of glucocorticoid (GC) therapy is GC-induced ocular hypertension (OHT) and GC-induced glaucoma (GIG). GC-induced OHT occurs in about 40% of the general population, while the other 60% are resistant. This study aims to determine the genes and pathways involved in differential GC responsiveness in the trabecular meshwork (TM). Using paired bovine eyes, one eye was perfusion-cultured with 100nM dexamethasone (DEX), while the fellow eye was used to establish a bovine TM (BTM) cell strain. Based on maximum IOP change in the perfused eye, the BTM cell strain was identified as a DEX-responder or non-responder strain. Three responder and three non-responder BTM cell strains were cultured, treated with 0.1% ethanol or 100nM DEX for 7 days. RNA and proteins were extracted for RNA sequencing (RNAseq), qPCR, and Western immunoblotting (WB), respectively. Data were analyzed using the human and bovine genome databases as well as Tophat2 software. Genes were grouped and compared using Student's t-test. We found that DEX induced fibronectin expression in responder BTM cells but not in non-responder cells using WB. RNAseq showed between 93 and 606 differentially expressed genes in different expression groups between responder and non-responder BTM cells. The data generated by RNAseq were validated using qPCR. Pathway analyses showed 35 pathways associated with differentially expressed genes. These genes and pathways may play important roles in GC-induced OHT and will help us to better understand differential ocular responsiveness to GCs.